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Abstract 

The master formula approach to chiral symmetry breaking proposed by Yamagishi and Zahed is 
extended to the Uij,(3) xUl(3) group in which effects of the Ua(1) anomaly and flavor symmetry 
breaking, ^ 7^ nig, are taken into account. New identities for the gluon topological 
susceptibility and tt^, rj, rj' — )• ^(*)-)/(*) decays are derived, which embody the consequences of broken 
chiral symmetry in QCD without relying on any unphysical limits. 

PACS numbers: ll.30.Rd, 11.40.-q 
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I. INTRODUCTION 



Zahed is a powerfu 

n 

doscalar mesons [5|, 



Phenomena of the rj and r]' mesons have been an attractive subject both theoretically and 
experimentally. Various properties of these mesons are closely related to the 11^(1) anomaly 
of QCD and the flavor mixing arising from the mass differences among u, d, s quarks. Since 
the Uyi(l) anomaly is a realization of the nontrivial topology of the gluonic configurations, 
investigating reactions and decaying processes associated with the r] and rj' mesons will be 
one of the most feasible approaches to accessing the fundamental aspects of QCD. 

Several theoretical approaches, such as the chiral perturbation theory [l, 2| and the zero- 
momentum Ward identities with PCAC hypothesis js, 4|, have been applied for analyzing 
reaction processes of t] and rj' mesons and exploring the gluonic context of the low energy 
QCD. Most of these investigations start from certain unphysical limits such as chiral, soft 
pion, and large Nc- 

The master formula approach to chiral symmetry breaking proposed by Yamagishi and 
l1 tool for analyzing hadronic processes which include ground state pseu- 
, |6| . The approach is based on a set of master equations, which fully takes 
account of the consequences of broken chiral symmetry without relying on any unphysical 
limits or expansion schemes. Also, the master equations provide a systematic procedure, 
called the chiral reduction formula (xR-F), to derive the chiral Ward identities satisfied by 
scattering amplitudes involving any number of pions with their physical masses. The ad- 
vantage of this approach is that one can investigate the reaction dynamics, which cannot 
be determined by chiral symmetry, separately from the general structure required by the 
broken chiral symmetry. Once such a separation is made, any models and/or expansion 
schemes can be employed for describing the reaction dynamics without any contradiction 
with the constraints from broken chiral symmetry. A number of investigations based on 
the approach has been carried out for hadron reactions in the resonance region Ulnl and 
hadronic matter 12l-ll8|. 

Up to now, the approach has been^ormulated within the two- flavor SU/j(2) x SUl(2) 
and three-flavor SUr(3) x SUl(3) [6] groups in the isospin symmetric limit. To analyze 
reaction processes involving t]' and address the features of the low energy QCD mentioned 
above, however, we need to extend the approach to the Uij(3) xUl(3) group incorporating 
the full flavor symmetry breaking due to m„ 7^ rrid 7^ and the U^(l) anomaly. In this 
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paper we will describe how to make such an extension. 

This paper is organized as follows. In Sec. |TT]we first review thejnaster formula approach 
for the isospin symmetric SUij(2) x SUl(2) group proposed in Ref. js] and outline a derivation 
of the master equations with this simplest case. Then in Sec. IIIII we describe the extension 
to the Uij(3) xUl(3) group including finite quark masses and the U^(l) anomaly. Several 
applications of the \J r{3)x'[J ^(3) master equations are presented in Sec. [IVl Summary and 
outlook are given in Sec. |Vl 



II. REVIEW OF ISOSPIN SYMMETRIC SU(2)xSU(2) CASE 
A. The Veltman-Bell equations 

Consider QCD with massive u and d quarks in the isospin symmetric limit: m„ = rrid = 
171. The Lagrangian can be written as 

J^Zb = [^qcb\o) + + l5a;)^q{x) - q{m + s - ip'^J^r^qix), (1) 

where ['Cqcq^](o) is the QCD Lagrangian in which the quark masses are set to zero; r° 
(a = 1, 2, 3) is the Pauli matrix for isospin; q{x) is the isodoublet quark field q = {u, d)^. The 
vector, axial- vector, scalar, and pseudoscalar external fields, (j) = (wj^, a^, are treated 

as sources to generate the corresponding currents and densities, O = (V^, A'^, E, 11"), which 
can be defined by O = 6{J d'^xC^^^)/6(f). 

A fundamental quantity in the theoretical framework developed in Ref. 5|] is the extended 
S-matrix operator S[(f)], a functional of the external fields (p. This operator is unitary, 
S'^S = SS'' = 1, and is related to the vacuum-to- vacuum transition amplitude in the presence 
of the external fields 0: Z[(j)] = (0 out|0 in)^ = (0 in|5[0]|O in) = (0 out|5[0]|O out). The 



Schwinger action princip 

o = {v;,aij:,u'^) as Is 



e allows one to express the quantum operators corresponding to 



and their T*-product as 



r.iofe)-o(.,.)i^(->r5'^...^5. (3) 



3 



A more detailed description of the theoretical formulation based on the extended S-matrix 



can be found in the literature 
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and this will not be discussed here. 
The system described by Eq. ([T]) and its effective theory have an approximate 
SUij(2)xSUL(2) chiral symmetry explicitly broken by the quark masses. It is known that 
such systems satisfy the divergence equations for the vect^ and axial currents, which, fol- 
lowing Ref. ^1, we call the Veltman-Bell (VB) equations 
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23| . The explicit forms are 



(4) 



V^^'^'A"^ + a^'^'V;^ - (m + s)W + p'^T. = 0, (5) 

where we have introduced the notation JC^^ = e°'^'^X^ (applicable to any quantity with 
one-isospin index X^); V"'^ = 5"^ (9^ + t;' 



ac 
■II ■ 



With Eq. ([2]) , the VB equations ^ and (|5]) can be rewritten as a set of linear equations 
of the extended S-matrix: 



5af, - dp"" 



yacj_ 



— [s + m) 



5 = 0, 
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6_ 

6^ 



S = 0. 



(6) 



(7) 



Applying functional derivatives of to Eqs. and ([7]) and using Eqs. ([2]) and ([3]), one can 
derive the vector and axial Ward identities satisfied by the operators O = (V^, A^, S, 11"). 



B. Master equations for the chiral symmetry breaking 

The VB equations (jl]) and ([5]) [or equivalently Eqs. ([6]) and (I7l)] are satisfied by the 
systems in both the Wigner and Nambu-Goldstone (NG) phases. However, if the chiral 
symmetry in the system is spontaneously broken, the NG bosons (pions) appear and couple 
to the axial current and the isovector-pseudoscalar density: 

(0|yi;:(x)|7r^(p)) = /.zp^5»^e-^P^ (8) 

(0|n'^(x)|7r''(p)) = G^^^e-'P". (9) 
Here f^^ is the pion decay constant, which remains finite in the chiral limit in contrast to 



those of the pion excitations 2J]; G is the pseudoscalar coupling constant. By taking the 



4 



matrix element of the axial VB equations ([5]) between the vacuum state (0| and the one-pion 
state |7r"(p)), and using Eqs. ([8]) and ([9]), we obtain the following mass relation: 

Uml = niG, (10) 

which shows clearly the well-known result that the nonzero quark masses are responsible for 
the nonzero pion masses. The master equations proposed in Ref. [5j] may be understood as 
the VB equations incorporating this information on chiral symmetry breaking. 

It was shown in Ref. [sl that the information can be incorporated into the VB equations 
by making the following modifications: 

1. Introduce new pseudoscalar and scalar external fields, J" and Y, defined by 

r = Gp- + UVl'a^^ (11) 

Y = Gs, (12) 
and treat = (a^, v'^, Y, J"-) as independent external fields. 

2. Introduce a new extended S- matrix S as 

S = Sexp{-i6I), (13) 

with 



SI = / d^x 



f2 

Y{x)G-^C + ^a'"'{x)a''^{x) 



(14) 



where a new constant C is introduced. (The physical meaning of C is explained below.) 

With these modifications, the new current and density operators, O = (jy^, J^^, <5", vr"), 
defined hj O = —i<S'^{6/6(j))<S, are related with the original current and density operators, 
O = (\/", A^^, S, n"), as follows: 

= ji, + fx - ^^r^^ (16) 

S = Ga + C, (17) 

n" = Gvr". (18) 

Here the new pseudoscalar density n"' satisfies {0\7c°'{x)\tt^{p)) = (5'^^e~*^^. This allows one 
to identify tt" with the normalized interpolating pion field. The change of the field variable 



jf- — )■ defined by Eq. (fTTj) is responsible for tlie separation of the one-pion component from 
the axial current and thus is introduced to indicate that the system under consideration 
includes pions. The new axial current has no one-pion component surviving on pion 
mass-shell, (0|j^^(x)|7r^(p)) = 0, in contrast to the original AJ^. 

The second modification defined by Eqs. f ll3p and f ll4p is introduced to indicate that 
the action describing the original S-matrix operator, iS, should contain contact terms, as 
defined in Eq. f lT^ . if the chiral symmetry is spontaneously broken. A justification of this 
statement has been provided in Ref. jsl making use of the gauged nonlinear sigma model. 
It is also noted that the difference between S and S just comes from the contact terms of 
the c-number external fields and thus does not affect the physical observables. 

The practical role of the second modification is twofold. The first term in the integrand of 
bX is introduced to take account of the quark- ant iquark condensation. In fact, it explicitly 
introduces a shift in the scalar density [see Eq. (fT7|l ]. amounting to the new constant, C. 
This constant carries part of the information on the condensation. (We allow (o") 7^ 0.) It 
turns out that C is expressed as the product of the pion decay constant and the pseudoscalar 
coupling constant: 

C = UG. (19) 

This follows from the fact that vr" is the normalized interpolating pion field. (See Appendix Rl 
for the derivation.) The mass relation ffTOl) can be rewritten as 

If {a) = 0, then C = (S) and Eq. reduces to the Gell-Mann-Oakes-Renner (GMOR) 
relation. Therefore, (a) represents the deviation of the mass relation from the GMOR 

n 

relation. In fact, the on-shell (V/tt) expansion scheme proposed in Ref. [5[, which is the 
expansion in l//,r around the physical pion mass and is constructed so that the GMOR 
relation holds at the leading order, leads to {a) = + C(/^^) and C = (S) + 0{f~^). 
The second term in the integrand of 6X results in the appearance of /^a^ in Eq. ( |T6|) . This 
ensures the existence of the contact term S"'^flg,i, in the two-point function of the axial 
current, 

z / d'xe^^-{T*[A';,{x)A''M]) = 5'^' fig,, - 5'^'fl^^ + • • • , (21) 

which must appear in the two-point function to have a correct chiral limit. (The symbol ( ) 
denotes the vacuum expectation value.) 
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Substituting Eqs. ([II]), (Il2]), and ([I5])-(II8]) into the VB equations and using Eqs. ([lO]) 
and f|l9|) . we have 

["^'""jW + ^'""fA, + r^'] = 0, (22) 



- J"(x) - K'iW'^'^f^^ + a^^^fy^) - [r - UV^''%]a. (23) 



These equations can be written in the functional derivative form as 



K"t- + <T- + J""tt- 1 ^5 = 0, 



SU(2);^j6 



5 = 0, 



with 



-"■SU(2) 



By introducing the retarded and advanced Green functions satisfying 



(24) 
(25) 



the axial VB equation ( 123|) can be formally solved for the interpolating pion field vr. The 
functional derivative form of the solution is written as 

5S 



^SnUx)+^S / (x, y)irsu(2)(2/)^.n(y) - / (x, y)i?^u(2)(2/)'5 



= z<(x)5 + ^ J AG^''(x,y)Ks%2)(y)vr(n(y)5-y AG^''(x,l/)/2^u(2)(l/)5.(26) 

Here ttj^ is the in-state asymptotic pion field, tt" -t- ttj^ + ■ ■ ■ (t — > — oo); -Rgu(2)(^) defined 
by 

' (27) 



R^lK9\ix) — -RcTTO-llx) + K, 



ab 



i'su(2)v-^; - -'i'su(2)v-^; -r -"'su(2)v-^;^jf,j^^^ 



with 



-^SU(2)(3^) 



1 A 



[X]. 



(28) 



Equations flM|) and fl26]) constitute the master equations for the SU/i;(2) xSUl(2) chiral 
symmetry breaking in the isospin symmetric limit. 
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C. Chiral reduction formula 



From the axial master equation (l26l) . we can derive the commutation relations between 
the creation and annihilation operators of the pion and the extended S-matrix S, 

[aUk),S] = i?gu(2)(fc)'5, [S,at^{k)] = RI^^,){-k)S, (29) 

where -Rsu(2)(^) ~ / d'^xe^'^^ R'^^^2)i^) ■ Note that we can rewrite the commutation relations 
in this form without using the asymptotic pion field. [Compare them with Eqs. (6.2) and 
(6.3) in Ref. 

Iterative use of Eq. (l29l) results in the fo^' the on-shell scattering amplitudes involving 
any number of pions with their physical masses, 

{a]kiai,- ■■ ,kmam\S\f3;hbi,- ■■ ,/„&„) 1^=0 = 

[-^su(2)(^i) ■ ■ ■ -^su(2)(^"^)-^su(2)(~^i) ■ ■ ■-Rsu(2)(~^")]s("l'^|/3)U=o, (30) 

where ki (/j) and Oj (foj) are, respectively, the four momentum and isospin indices of the 
outgoing (incoming) pions; a and P stand for states of other particles. Here we consider 
the case that no two pions have equal momenta. The symbol [ ]s represents normalized 
symmetric permutations of the functional derivative operators, 

[V,---Vr,]s = ^_ 1^1- --l^n. (31) 

perms. 

This operation shows clearly the crossing symmetry in Eq. (!30|) . By using Eq. ( 130|) to- 
gether with Eq. (|3]), scattering amplitudes are expressed in terms of Green's functions of 
the operators O = (jy^, j^^, o", vr"). The xR-F takes the form of functional derivatives, and 
all constraints which stem from broken chiral symmetry are contained in Rsuiaxi^)- The 
extension of the xRF to the off-shell pions has been discussed in detail in Ref. |25|. 

III. EXTENDING TO U/j(3)xUl(3) WITH FLAVOR SYMMETRY BREAKING 

In this section, we extend the master equations reviewed in Sec. [Ill to the U/j(3) x 11^(3) 
group with rriu 7^ 7^ 171^. In the remainder of this paper, the term "pion" ("vr") is used for 
expressing the nonet ground pseudoscalar mesons generically, and the symbols 7r^'°, K'^''^, 
K~'^, rj, and rj' are used for referring to the specific mesons. 
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A. The U/j(3)xUl(3) VB equations 



The Lagrangian now includes massive m, d, s quarks: 

^QC?i = [^qcd]{0) + + l.a^^q - q{m-^ + - ip^^,)\^q - Ou. (32) 

Here = {u, d, s); m^A" = diag(m^i, m^, m^); 6 and u are the vacuum angle and the gluon 
topological charge density, respectively. The term [>Cq[?^](o) represents the QCD Lagrangian 
with the quark masses, the vacuum angle, and all external fields set to zero. The fiavor 



matrix A" is taken to be one of the Gell-Mann matrices for a = 1, ■ ■ ■ , 8 and A° = a/2/31 so 
that they satisfy Tr[A"A^] = 25"''. The vector, axial- vector, scalar, and pseudoscalar external 
fields and the vacuum angle, = (f^, a^, s°,p", — 6*), are sources to generate currents and 
densities, O = (V^", A^, S*^, E", cu), which are obtained from O = 6{J d'^xC^^j^) / dcj). As in 
the SU(2) case, the quantum operators corresponding to O are defined by Eq. ([2]). 
The VB equations for the Ui(;(3) xUi(3) group can then be written as 

V^^^'A"^ + ^llVl + ^^Y}' - (s + mg)""n" - Tr[A^u; = (34) 

Here we have introduced the notation^ X"*^ = d°^'^X^ and X""" = f"-^^X^ (applicable to 
any quantity with one-fiavor index X*); V^'^ = dp^d""^ + w^'^. It is well known that an ad- 
ditional nonconserving contribution in the axial VB equation, which amounts to — Tr[A'^]a;, 
is attributable to the Ua(1) anomaly. Furthermore, for the Uij(3) xUl(3) group, the non- 
Abelian anomaly fi" associated with the external gauge fields and aJJ also appears. Its 
explicit form is |26| 



= — ^^^'^''"Tr 



2 \ fJ-i^ pa ' 2 /.w pa 



8 32 



(35) 



with = dpVy - dyVp - i[vp,Vy\ - i[ap,ay\, = d^ay - d^a^ - i[vp,ay] - i[ap,Vu], 
Vp = f"(A"/2), and = a'^(A"/2). We take a convention with £0123 = +1- 



1 The structure constants and d"'"' are defined by = -(i/4)Tr[[A'', A'']A'=] and d"^^ 

(l/4)Tr[{A'^,A''}A=], respectively. Witli tiiis definition, we have f"" = and d°'"' = ^2735''^ 
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The functional derivative form of the VB equations is given by 

5 = 0, (36) 



B. Flavor symmetry breaking and U^(l) anomaly 



S = iQ^'S. (37) 



Besides their key role in the singlet sector, the mass differences among the u, d, s quarks 
and the U^(l) anomaly are responsible for the major complication in the construction of 
the U f({3)x\J 1^(3) master equations. Therefore we first examine those effects carefully. 

The mass differences of the quarks violate the flavor symmetry explicitly as well as the 
chiral symmetry. Thus the flavor U(3) basis cannot be the mass eigenstates of hadrons and 
mixing of the U(3) basis occurs. Necessary information on the flavor symmetry breaking 
here is summarized in the following matrix elements, 

{0\A';,{x)\P{p))=fftp,e-'^^, (38) 

(0|n"(x)|P(p)) = G^^e-'P"", (39) 

{0\u{x)\P{p)) =A^e-'P'', (40) 

where P denotes the physical pion states (the mass eigenstates): P = (vr^'", K^'^, K^, t], rj'). 
(Capital indices are used to represent the mass eigenstates.) Because of the flavor symmetry 
breaking, the pion decay constant /^^ and the pseudoscalar coupling constant G""^ now have 
two indices that specify the flavor U(3) basis and the mass eigenstates^. Also, it is noted 
that has a nonzero value only for P = n^, rj, rj'. 

From the axial VB equation ( !34l) and Eqs. ( l38|) -(l40l). we can derive the U(3) version of 
the mass relation: 

ff^rnlfP = rrf^'G'"' + Tr[Al^^, (41) 



The breaking pattern of the flavor symmetry, i.e., m'^X"- = rn^A" + m^A^ + m^X^, imphes that and 
G"^ are block diagonal in (1, 2) x (7r+, tt"), (4, 5) x i^-), (6,7) x {K° , K°), and (0, 3, 8) x (tt^, 77, r/'). 
In the isospin symmetric limit they become (1,2,3) x (7r+, tt^ , 7r°), (4,5,6,7) x {K'^ , , K'^ , K*^), and 
(0,8) X (77,77') and the first two blocks can be diagonalized. 
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where is the diagonal mass-squared matrix of the physical nonet pions. The appearance 
of is attributable to u in the axial VB equation and thus is nothing but a consequence 
of the U^(l) anomaly. 

As already mentioned above, the U^(l) anomaly leads to the nonconserving term pro- 
portional to u in the singlet axial VB equation. Because of this, the operators A^^ and 



Lo get renormalized under a change of the QCD scale [l]: (A°)ren. = ^aA^ 



LO + {Za — l)(l/Tr[A''])(9^74°, where Za is the renormalization factor. Accordingly, the 
singlet pion decay constant and A^ also get renormalized as 

(/Dren. = ^a/^, (42) 
(i^ren. = + {Z A - 1) ;^/°^(m^)«^. (43) 

These are thus scale dependent and not physical constants. This is in contrast to the octet 
pion decay constants, /"^, with a 7^ 0, which are scale independent. It is noted, however, 
that the combination f^'°^{m'^)^^ — Ti[\^]A^ is invariant under a change of the QCD scale. 

As pointed out by Shore and Veneziano (see, e.g., Refs. {4, 27, 28|), the renormalization 
group (RG) variant /°'^, which is defined as a coupling of the singlet axial current to the 
physical pions, does not satisfy GMOR type mass relations. They introduced a RG-invariant 
decay constant f^^ and showed that the use of this new constant provides a natural extension 
of the GMOR relation to the singlet sector. 

Following Shore and Veneziano, let us introduce the RG invariant /°^. Defining 5/°^ = 
f^^ — f^^, the combination {fT^m\Y^ — Tr[A'']^'^ can be rewritten as 

CUmir - Tr[A°]^^ = (/.m^)°^ - Tr[A°][^^ - (^/.m^)^- (44) 

Because (/tt^^)"^ — Tr[A°]yl^ and (/^m^)°^ are individually RG invariant, A^ = [A^ — 
{6 fT^ml)^^] is also RG invariant. The values of and A^ must be extracted from the 
experimental data. With this modification, the mass relation ( I^Tl) becomes 

iUrnlfir = {m,Gfjr + TT[X'^]{AfJf 

= {m,Gf:t + (45) 



where ff = ff for a 7^ 0; (/J) = {f^Y" is the transpose of ff- A"^ = Tr[A"](^/^ 



T\b 
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C. Master equations for the Uij(3) xUl(3) symmetry breaking 



We can incorporate the information on the symmetry breaking described in Sec. IIIIBI in 
a completely parallel way to the SU(2) case: 

1. Introduce a new pseudoscalar and scalar external fields and defined by 

yP = gaQaP^ (^7^ 

with al = al + 5"°(l/Tr[A°])9^^, and then consider cf) = {a"^, v^, Y^, J^, -6) as inde- 
pendent external field variables. 

2. Introduce new extended S-matrix by 

5 = 5 exp [-i {61 + 6In)] , (48) 

where 

SX = I d'x {yP{G-')P^C^ + \a^''{U^r~al 

5Xn = ^.e^'""' I d'x (^a;vrapV^S° - ^a°9.aJv^S°^ . (50) 

Note that the new constant C"^, which corresponds to C in the SU(2) case, now has a 
fiavor index. 

With these modifications, the new current and density operators O = (jy^, j^^, o-'^, tt^, a)) 
are defined hj O = —iS^S/Scp with (p = {v'^,a1j^,Y^,J^,—9). The relations between the 
new and original operators are given by 

= Jy^ + (/.vr) + -j-j^, (51) 



and 



(52) 

ir[/vj uw — 

_ G^'^n^, (54) 
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+^^Wjy''Tal + + (Sfja"'n]e + (55) 

As in the SU(2) case, we can identify the new pseudoscalar density vr^ with the (normahzed) 
interpolating pion field satisfying (0|^^(x)|P(A;)) = 6^^e~^^^. 

The modifications induced by Eqs. (H6|) . (1471) . and fHOj) can be understood as a straight- 
forward extension of, respectively, Eqs. flTTl) . f[T^ . and f[T^ in the SU(2) case. With the 
modification (l46l) . the one-pion component is separated from AJJ and u, and the resulting 
new operators and have no one-pion component surviving on pion mass-shell, i.e., 
(0|j^^|P(p)) = {0\u\P{p)) = at p2 = {ml)^P. From Eq. ([53]), the new constant 
again carries part of the quark-antiquark condensate: (S") = G"'^{cr^) + C". The terms in 
Eq. poj) . except for the first one, are introduced such that the two-point functions including 

and/or u have the correct chiral limit in the presence of the pions. 

A comment on the term 6Xq is needed since it does not appear in the SU(2) case. As 
discussed by Kaiser and Leutwyler jl]], the non-Abelian anomaly, fi'*, which includes the 
external singlet axial gauge field a°, is not invariant under a change of the QCD scale 
because of the Ua(1) anomaly. This leads to an inconsistency with the RG invariance of the 
VB equations^. The term SX^ is introduced to cure the inconsistency and corresponds to 
the sum of the two contact terms. Pi and P2, of Eq. (78) in Ref. With this additional 
term the non-Abelian anomaly fi" is replaced with the RG invariant Qq, which is of the 
same form as fl"-, but a° is replaced with — (l/Tr[A°])9^6'. 

Taking account of the modifications described above, the VB equations fl33l) and flMl) can 
be rewritten as 

V-^'fy^ + + {G-'jf iGn)" + [ {G-'Yf + rrz^^"] {Gaf + Xv,^'' + Xv2 = 0, (56) 

and 

-ab 



iL-'riG-'J) iGaf + {f-Y^{Q^,-xl). (57) 



^ The RG invariance of the VB equations arises from that of the S-matrix S. 
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Here we have introduced 



= atCUEt~^"' - afiUfjr^^d^e + {G-'Yf c\ (59) 



^ab ^ _ (^gl^a^yb + {G-^Y) {Gf;^^ - (61) 

J'' = J""- ifif^'i^xr + [(/:)™° + if^y''^''']Tf^f- (62) 

The terms xvi, Xv2, and xa are attributable to the flavor symmetry breaking and thus they 
vanish in the flavor symmetric limit. 

In deriving Eqs. ( 156|) and ( 157|) . we have made use of a relation satisfied by the new 
constant C", 

C"' = {GfJr\ (63) 

which can be regarded as the U(3) version of Eq. (fT9|) . This relation follows from the fact 
that TT^ is the normalized interpolating pion field and can be derived by the same strategy as 
described for the SU(2) case in Appendix |Al Then the mass relation (H5|) can be rewritten 
as 

iUmlf^r" = {m,Cr + A^". (64) 

From {mgC)^ = {rrigC) and A'''' = for a ^ [recall that A''^ = Tr[A''](^/J)^], one can 
fix A""^ up to one constant: A"'^ = Tr[A"]Tr[A^]y4^. Within the l/fn expansion scheme jsl, 
we have (a^) = at the leading order, and therefore C° = — (a/2/3)(('um) + {dd) + {ss)), 



G^ = -{{uu) - (dd)), G^ = -{^/T/3){{uu) + (dd) - 2{ss)), and G" = for the other flavor 
indices. In this case Eq. ( IMl) reduces to the generalized GMOR relation given by Shore jj], 
which is derived by making use of the Ward identities in the zero-momentum (soft pion) 
limit. Because the leading order of the 1/fn expansion scheme corresponds to taking the soft 
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pion limit, we can identify the constant with the nonperturbative coefficient appearing 
in the gluon topological susceptibility in QCD j4|. The quantity (a^) again represents the 
deviation of the mass relation ( I64p from the GMOR relation in the same way as for the 
SU(2) case. 

The functional derivative form of Eqs. (156 p and fl57p is given by 



mx)s = 0, 

n(x)s = 0, 



(65) 
(66) 



with 



+ <(^)Yzr7z:^ + (G-'Jr\x)G 



bP 



5 



(5<(x) 



+ [(G'-^r)"'(x) + mf]G'^ 



+ xfAx 



5Yi 



\x 



6JP{x 
6 



6JP{x) 



+ iXv2ix), (67) 



Tlix) = ff [-n5^^ - K)^«] jj^^ + ffR'ix 



(68) 



6JQ{x] 



+ iJ^(x) 



+(/; 



l\Pa 



6 



+ a 



yab " I ab_ 

^ 5a^^{x)^->^ 5vl{x) 



+ Tr[A''] 



8e{x) 



+{f;Y''{G-'jf\x)G'^ 



6 



Kf;;'rmx) - xa{x)], 



6YQ{x] 

By introducing retarded and advanced Green's functions satisfying 

[-nj^^ - (mlfQ - {L-Y^K^\x)ft^] GZ{x,y) = 6''^6\x - y), 



(69) 



one can formally solve Eq. (I66l) as 

6 



5JP{x) 



S = ^S7rax)+^S I d'yC^J^ {x,y){f-Y''K^'iy)fn''<{y) 



l\Qa T^ab/„\ fbRRt 



-J d''yG'^'^{x,y)R'^{y)S 

<(x)5 + ^ J d'yG'^'^ix,y)ifYf'^K-\y)C7r^M^ 
d%G^^'^ix,y)R'^iy)S, 



(70) 
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with 



6 



(71) 



6JQ{xy 



where is the in-state asymptotic pion field. Equations ( l65ll and (!70l) are the desired 
extension of the master equations to Ui?(3) xUl(3), incorporating the finite u,d,s quark 
masses and the Ua(1) anomaly. The major consequences of fiavor symmetry breaking are 
the two- index character of the various constants and the appearance of Xvi'^^^ Xv2i ^^"^ X'a 
terms. On the other hand, those arising from the 11^(1) anomaly are the mass relation ( 16^ 
and the appearance of the operator Cj. 

D. Chiral reduction formula 

The commutation relations of the pion creation and annihilation operator with the ex- 
tended S-matrix S are given by 



where R^{k) = J dxe^^^ {x) . The Ur{3)xUi[3) version of the xR-F fo^' on-shell pions can 
be expressed as 



where fcj (/j) is the four momentum of the outgoing (incoming) pion Pj (Qi), and a and f3 
label states of other particles. Here we again consider the case that no two pions have equal 
momenta. The symbol [ ]s is defined in Eq. ( 13T|) . 

Before closing this section, we note that the singlet axial current j^^ (or its functional 
derivative form) always appears as the RG invariant combination d^j\^ — Tr[A°](u) in the 
axial master equation flTOj) and the xR-F fTTSl) . [In the absence of the external fields, j^^ 
and w° are renormalized as (j°^)rcn. = ^aJa^ and {Q)ren. = w + {Za - l)(l/Tr[A°])9^j^^, 
respectively.] Therefore, the existence of w, which originates from the U(l)yi anomaly, is 
crucial for ensuring the RG invariance of the master equations and the xR-F- 



KiklS] = R''{k)S, [S,a^J{k)] = R^{-k)S 



(72) 



{a■Pl{k^),■■■ 



Pm{km)\S\f3;Qi{k),--- ,Qn{ln))U=0 = 

[P^H^i) ■ ■■R''-ikm)R'^'i-k) ■ ■ ■i?^"(-/n)]s(«|5|/3)U=o, (73) 
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IV. APPLICATIONS 



As an illustration, we will present several applications of the U/j(3) xUl(3) master equa- 
tions and the x^F- 

A. Gluon topological susceptibility 

By using Eqs. (E]) and (133]) . we can derive the chiral Ward identity for the gluon topological 
susceptibility x'- 



X = -ij d'x{T*[uj{x)ujm) 
= A^-i j c/^x(T*[w(x)w(0)]) 

-t2VQA^{f;Y' I d^x{T*[coix)n^m)- (74) 

This shows clearly how the pion poles contribute to the gluon topological susceptibility. The 
constant A^ of the first term corresponds to the leading contribution in the large Nc limit 
of QCD with massive quarks, x = A^ + 0{1/Nc) |28|, The appearance of A^ is ensured 
by the modification introduced by the third term of 6X [Eq. fH9|) ]. The RG transformation 
property of u implies the RG invariance of the zero-momentum projected two-point functions 
/ d'^x{T*[Ld{x)fc^{0)]) and / d'^x{T*[uj{x)uj{0)]) , and therefore Eq. ([74]). 

Making use of the chiral Ward identities of (T*[a)(x)7r^(0)]) and {T*[7r^ (x)?!-^ (0)]) derived 
from the axial master equation fITOl) (see Appendix[B]for the results), the above identity can 
be further written as 

X = A,-QAl[{UmlfJ)-T 

-^{l-QA^[{^mlfJ)-Tr I d'x{T*[u;{xMO)]) 

■cb 



+QA'A(f^mifi) rH^{G{a)) [(/.m^/j )"?"• (75) 

Here the second term comes from the pion pole in (T*[7r^(a;)7r'^(0)]). This is the most general 
expression constrained only by the broken chiral symmetry. 

We observe that our result consistently reduce to those obtained in previous works by 
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taking appropriate limits. With the mass relation (IMj) . we obtain 



\ q=u,a,s ^ ^ / 



-1 



(76) 



where Cu = C + C3, Cd = C - C3, and = (Co - V2Cs)/VQ with C = (2Co + Cs)/V^. 
The on-shell expansion scheme [s] gives Cq = ~{qq) {q = u,d,s) at the leading order. In 
this case, the first and second terms of our general expression (175|l (i.e., the contributions 
from the leading term in the large limit and the pion pole term) reproduce the classic 
result ^, l29|. Also, our result approaches zero in the chiral limit, x {rriq — )■ 0), showing 



the noncommutative character of the Nc ^ 00 and ruq — )■ limits 



30|. 



The third and forth terms in Eq. f l75p are new. The third term shows how the higher 
meson states X, including hybrids states and glueballs, contribute to the gluon topolog- 
ical susceptibility through (0|a;|X) 7^ 0. The forth term is proportional to (a^), which 
characterizes the deviation of the mass relation (EH) from the GMOR relation. 



B. Two-photon decay of vr*^, 17,17' mesons 



Next let us consider the two-photon decay of the 7r°, rj, and t]' mesons: P{p) — )■ 
'y^*\qi)l^*\q2) with P = {n^,ri,ri'). For this purpose, we need to evaluate the amplitude 
J dXq^dyq^{0\ST*[j^^{x)j^^{y)]\n^{p))\4,=o, where J dxq = J d^x exp{iqx). The xRF gives 

dXqMO\^T*[fy^{x)j'y^{y)]\n''{p))U=o 



dxq^ dyq2 \ 



6 



6 



5<(x) 6v^^{y) 



i?^(p)(5)| 



(j)=0 



dxg.dyq^dz^pil 
where = d^j 



-l\Pc 



6 



6 



_ 6v^{x) 5vl{y) 



(77) 



Tt[\^\Cj. With expressions for the electromagnetic current, 

,8 



em 



UvY + {^/^){jvYi ^iid the external electromagnetic field, —eA^ = v 



we can 



derive a general expression for the P — t- decay amplitude: 

9-np{ql^ql]P^) = 



l\Pc 



em g^2 



with (Cgjj^, Cgj^, Cgj^ ) = (2/3,2^3/9,476/9). Here we have defined 

dXq,dyqM^T*[JUx)J:^{y)]\'r^''{p))U=o = ^e^''''^qlaq2^ 



(79) 
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(80) 



and 



dx,,dy,,Jz_,{T*[ji:jx)j:Jy)dxjfiz)]) = 

6^^"^i„g2/.[F^Vi(^i'^2;p') - F^,A2i(llQhP% (81) 
modulo {2tt)'^6^{p — qi — q2)- In Eg- (1ST]) we have used the general expression of the vector- 



vector-axial correlation function 



31 



with 



32| [denoting .-{qlqlp^) as F: 



] I dx,,dy,,{T*[jUx)j:jy)jfm) = i2^yy^^^^^ 



i=l 



e^'""'qia, 



-'2 



Xfia/B V rXfiu 

'?la'?2/35'2? -'4 



-'5 — e '?la'?2/35'l ; -'e 



Xfii/a 
e q2a, 



Xuaf} ^l 
e '?la92/35'l , 

e qiaq2i3q2- 



The broken chiral symmetry relates the P — )■ 7(*)7(*) decay amplitude to the correlation 
functions of the current and density operators ( !80l) and (ISTj) . The correlation functions 
themselves, however, cannot be fixed from symmetry requirements and thus some dynamical 
inputs are needed to evaluate them. 

Equation fl78|) is applicable both for on- and off-shell pions. Using the PCAC hypothesis 
5'77p(Q'i5 '?2j "^p) ~ Q'yypilhlh^) Setting qi and q2 to the photon point, qf = q^ = 0, 
our general expression (178]) consistently reduces to the two-photon decay formula of 7r°, rj, rj' 
mesons derived by Shore j^, |28[ , 



l\Pc 



87r2 



V6F^^<i(0,0;0) 



(82) 



if we identify F^^(2,(0, 0; 0) as Agc-y^- Here we have used P^^^i (0,0 0) = F^V2(0'0;o) 







32| . The RG invariance of the operator = Tr[A'^](u) implies that the combination 

'\/6p77(i — P°^yii + P°^A2 is also RG invariant. Therefore F^^i(0, 0; 0) is by itself RG invariant. 



V. SUMMARY AND OUTLOOK 

We have derived an extension of the master equations for chiral symmetry breaking 
proposed in Ref. to the Uij(3) xUi(3) chiral group, carefully taking into account 
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the Uyi(l) anomaly and full flavor symmetry breaking ruu 7^ rud 7^ m^. With the master 
equations and the xR-F; iiew chiral Ward identities for the gluon topological susceptibility 
X and P — )■ decay amplitude have been derived, showing how the constraints from 

broken chiral symmetry enter into those quantities without relying on any unphysical limits. 
Then we have seen that our general results consistently reduce to those obtained in previous 
studies by taking appropriate limits. 

The xR-F is applicable to any reaction processes which include ground state pseudoscalar 
mesons, e.g., vr, K, K, 77, 77' production reactions on a baryon target and heavy meson decays 
such as J/^ — 7- 3P, 72P, 73P, ■ ■ ■ with P = tt,K, K, rj, rj'. The heavy meson decays are inter- 



esting in relation to new meson resonance states appearing in the decay processes [33|, |34 |. 
A careful treatment of the flnal state interactions in the decay processes will be vital for ex- 
ploring properties of such new meson states. The enables one to separate details of each 
reaction mechanism from the general framework required by broken chiral symmetry, and 
thus will provide a useful theoretical basis for the analysis of such processes. Investigations 
in this direction will be discussed elsewhere. 
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Appendix A: Derivation of Eq. (I19p 

In this Appendix we describe the derivation of Eq. f|T9|) . which relate the new constant C 
to the pion decay constant and the pseudoscalar coupling constant G. The strategy used 
here can be straightforwardly applied to the U(3) case [Eq. ( 163|) ]: in this case the algebraic 
operations just become more complicated. 
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Originally the operator -Rsu(2)(^) following form 

-^SU(2)(^) = 



'SU(2) 



/ T JO, < ^ +0 ]^o.b ^ 
— J +-t^ ^SU{2)^ 



' fJSY 



+i 1 



GC 



[X . 



(Al) 



Equation ( 1191) is obtained by making use of the fact that the ti"" is the normalized interpo- 
lating pion field satisfying (0|'7r"(x)|7r^(p)) = §"-^e-^P^_ By applying the xR-F to the matrix 
element (0|7r"(x)|7r^(j9)), we have (note that S\0) = S^O) = |0) following from the stability 
of the vacuum state), 



(A2) 



with 



-p)S = + A^'^f'e-^P- + ^^afeg-ipx 



+ j- I d%e-^PyitprT*[fc'^ix)f^^iy)] + Oi<l>) 



d'ye-'^p^3\(y)7:l{x) 



+ j^e^''j d'ye-^^yAn{x - y)[-2^p^ + {dyr]fy,{y) 

-1 J d'yAnix-y)idyrsnsrA,iy),a^M + OmA3) 

where {dyY = [d/{dyi,)]. In the last step of Eq. (jMl), we have used Eqs. ([26]) and (ETj) with 
Rsu{2) replaced by Eq. (lAip . and the relation 

. 6 



—I 



Jtt 

-jf I d'ze-^^^T*[j%^{y)jU^)] + Oi^). 
Noticing that (OIj^^Jtt^) = 0, (0|j{^^|0) = 0, and = ml, Eq. (jM]) gives 

-| fir 

~ GC' 



(A4) 



(A5) 
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Appendix B: Commutation relations and chiral Ward identities 



1. Commutation relations 



The functional derivative operators defined in Eqs. fp7|) and fIBSl) satisfy the fol- 
lowing commutation relations [defining Ty{k) = J d'^x exp{ikx)Ty{x) and T^{k) = 
J d'^x ex]){ikx)T%{x)\. 

[T^{k),T^{k')] = -r^^T^ik + k'), (Bl) 

[T^{k),T'^{k')] = -r'^T^ik + k'), (B2) 

[n{k),T'^{k')] = -r'^T^ik + k'). (B3) 
With Eqs. ^ and dBS]) we further obtain 



[n{k),R^{k')] = -i^f-^^ r^^T^i^k + k' 
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iff + d""^^- 



6 



(B4) 



6YC{x) 

where d"'^'" = (i°^^(G'~^)^*G^*". If k'^ is the on-shell pion momentum satisfying (A;')^ = 
(m^)^^, we have \T^(k), R^{k')]S = 0. Therefore, we can make use of the same prescription 
as proposed in Ref. [25| to derive the off-shell extension of the U(3) (C3I)- 



2. Chiral Ward identities for the two-point functions 



Here we summarize several chiral Ward identities used for obtaining the results in 
Sec. HVAl 

Making use of Eq. ( l55l) and Eqs. ([2]) and ([3]), the two-point function of u can be expressed 

as 



J d''xe-'P^''-y\T*[uj{x)uj{y)]) = j d 



{T*[uj{x)u{y)]) 



+ i 



+ Ar^(p') j d''xe-'P'^^-y\T*[u{x)7c''{y)]) 
+ A/'^'b') J d^xe-'P^''-y\T*[oo{y)n''{x)]) 
+ A/'^(p')A/'^(/) [ d^xe^'P^''-y\T*[Tr''{x)n^{y)]),{B5) 
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with M^{p^) = y6A^(/-^)™ + {p^/VQ){Sf^f^. This identity reduces to that of the gluon 
topological susceptibility [Eq. ( !74|) ] by setting = and y = 0. 

Similarly, with the axial master equation f l70|) . the two-point function of the interpolating 
pion field can be written as 



+1 j d^xe'^^^-y^Ajf{x-y){f^^f'd'''''{a'') 
- j d'xe'^^^-y^ I dzAjfix - z)if;Y^T*[W^iz)7,''iy)]) 
= i j d'^xe'P^'^-y^Aj.^ix - y) 
+ i j d^xe'P^'=-y^Af^\x-y){f;Y"d''^^{^^) 

+ I J d^xe'P^^'-y^ j dzAjf{x~z)Al^{y-z) 

+ j rfW^^^-^) j dzdz'Af{x-z)Al'^{y-z') 

x{f~'fV^Y\T*[W\z)W\z')]). (B6) 

Here we have used W"" = d^j'X^ — Tr[A"]a); A^-^(a; — y) and A;^^(x — y) are respectively the 
retarded and Feynman propagators satisfying [—U^S^^ — (m^)^^]A;^^^(x — = 5^^5^{x — 
y). In the last step, we have made use of Eq. f lTOj) and the following Ward identities: 



d'xe^^^^-y\T*[3\^{x)ix\y)]) = j d'xe'^^^-y\3\^{x)nt{y)) 

+ j ciW^'(^-^)p^A^^(y-x)(/~i)^^(//)^"c/^^^(a^) 

d^^e^Pi^-y) [ dzA^^{y-z){f;Y''{T*[rA,{x)W'^{z)]), 



(B7) 
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- I 



V6 



rf4^e^p(x-.) / dzAf{y-z){f~^f\T^\Gj{x)W^{z)\). 

(B8) 



Setting p = and y = 0, Eqs. (IB6P and (IBSP become 



+^[(/X/J)-Tl(/.</J)-T^ 
+6[(/X/J)-T1(/.-^/J)-T° 

X /" c/^x(T*[w(x)w(0)]), (B9) 



and 



-VeliUmlfJrT J d'x{T*[u;ixMO)]), (BIO) 
respectively. Substituting these equations into Eq. ( [7i|) . we obtain Eq. ( !75|l . 
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